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Abstract 



We derive a new representation for spin-spin correlation functions of the finite 

XX Z spin-1/2 Heisenberg chain in terms of a single multiple integral, that we 

^+ ' call the master equation. Evaluation of this master equation gives rise on the 

^D • one hand to the previously obtained multiple integral formulas for the spin-spin 

correlation functions and on the other hand to their expansion in terms of the 

form factors of the local spin operators. Hence, it provides a direct analytic link 

^~H. between these two representations of the correlation functions and a complete 

(-H I re-summation of the corresponding series. The master equation method also 

allows one to obtain multiple integral representations for dynamical correlation 

functions. 
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1 Introduction 

One of the central question in the theory of quantum integrable models ^ |21 El 13 i^ ^^^ 
exact computation of their correlation functions. Apart from few cases, like free fermions 
[3 El 13 13 El ^1 or conformal field theories ^J, this problem is still far from its complete 
solution. In particular, the computation of manageable expressions for two point functions of 
local operators and their asymptotic behavior at large distance is a central open problem. If 
one considers the case T = 0, such a problem reduces to the calculation of the average value in 
the ground state \uj) of the product of two local operator 9i,92'- 

912 = (a;|^1^2|w). (1.1) 

There are basically two main strategies to evaluate such a function: 

(i) to compute the action of local operators on the ground state ^1^21'^) = P) and then to 
calculate the resulting scalar product gi2 = {oj\Cj); 

(a) to insert a sum over a complete set of states \uJi) (for instance, a complete set of eigen- 
states of the Hamiltonian) between the local operators 9i and 62 and to obtain the representation 
for the correlation function as a sum over one-point matrix elements (form factor type expansion 

iniiiniiiii) 

gi2 = ^^{id\9i\iJi) ■ {uii\e2\i^)- (1.2) 

i 

The aim of this paper is to give a direct, analytic relation between the approaches (i) and 
(ii) in the case of the XX Z spin- 2 finite chain. For the sake of simplicity we mainly consider 
the case of zero-temperature and equal-time a^ correlation function. It will be clear however 
that our method is more general. In particular, it allows one to compute dynamical correlation 
functions. The corresponding results will be described in a sequel to the present paper. 

We consider the periodical XX Z spin- 2 Heisenberg chain in an external magnetic field. The 
Hamiltonian of this model is given by [TS] 

M 
H =^ (cr.^cr^+1 + o-|^cr^+i + A(o-^o-^_^i - 1)) - hS-^, (1.3) 

m=l 

where 

1 ^^ 
Sz = i^Y.<^ \H,S,\=^. (1.4) 

m=l 

Here A is the anisotropy parameter, h the external classical magnetic field, a^'^ are the spin 
operators (in the spin- 2 representation), associated with each site of the chain m, and M is even. 
The quantum space of states is 7^ = 'S>^^iHrn, where Hm ~ C^ is called the local quantum 
space at site m. The operators Um^'^ act as the corresponding Pauli matrices in the space Hm 
and as the identity operator elsewhere. 



The method to compute eigenstates and energy levels (Bethe ansatz) of the Hamiltonian 
^l.'A\f was proposed by Bethe in 1931 in J^ and developed later in |17|ll8| [Tn]. The algebraic 
version of the Bethe ansatz was created in the framework of the Quantum Inverse Scattering 
Method by L.D. Faddeev and his school j2Ul 1211 Hj. However the knowledge of the correlation 
functions of the XX Z chain has for a long time been restricted to the free fermion point A = 0. 

In the case of the XXZ chain the approach (i) leads to multiple integral representations 
for the correlation functions. In the thermodynamic limit, at zero temperature and for zero 
magnetic field, such representations were obtained from the g-vertex operator approach (also 
using corner transfer matrix technique) in the massive regime A > 1 in 1992 J22j and conjectured 
in 1996 j23j for the massless regime — 1 < A < 1 (see also 24 ). A proof of these results together 
with their extension to non-zero magnetic field was given in 1999 j25| I26j for both regimes using 
algebraic Bethe ansatz and the actual resolution of the so-called quantum inverse scattering 
problem |25[ I27j . Using these results, the spontaneous magnetization was obtained in j2Sj . For 
the case of the XXX model (A = 1) this type of representation was later studied in j29l l3L)j . 
Integral representations for spin-spin correlation functions were derived in |311 I32j . Recently 
the generalization of this result for finite temperature was given in j^ . 

In the framework of the approach (ii) integral representations for the form factors of the 
XXZ chain in the thermodynamic limit were obtained in |221 1241 1^ 1351 OHl I37j . Determinant 
representations for the form factors of the finite chain were computed in |25[ 128j . An effective 
summation of the form factor series in the case of free fermions was done in |38[ 139j . 

To explain our method we will study the two point correlation function of the third com- 
ponents of spin (afa^^i). Following the papers j^ OHl 1^ we use for this purpose a special 
generating function {Ql^) (see (|4.1|) . (|4.2|) ). On the way to relate the previously obtained 
multiple integral representation of the spin-spin correlation function to its form factor type ex- 
pansion we derive what we call the master equation (see ()4.4|) ) for this generating function. It 
is given as a single multiple integral of Cauchy type over a certain contour F in C^ 

(Qlm) = i dzi... dzN F,{{z}). (1.5) 

r 

The integrand -Fk({z}) is a function of the N variables zi, . . . , zn- It is mainly given in terms 
of the eigenvalues r^ of the twisted transfer matrix (see (|3.H) ) and of the functions y^ defining 
the corresponding twisted Bethe equations (see (|3.5|) . (|3.()j) ). It is a periodical function of each 
argument Zj, vanishing at Zj -^ itcxD. Therefore there are two ways to evaluate the integral 
p.5() : either to compute the residues in the poles inside F, or to compute the residues in the 
poles within strips of the width iir outside F. 

The first way leads to a representation of the correlation function (cfc^+i) in terms of the 
previously obtained j31j tti- multiple integrals. Evaluation of the integral (|1.5|) in terms of the 
poles outside F gives us the form factor type expansion of the correlation function (i.e. an 
expansion in terms of matrix elements of a'^ between the ground state and all excited states). 



We would like to stress that the objects entering the above master equation for the finite 
XXZ chain are quite generic in the context of quantum integrable models solvable by the 
algebraic Bethe ansatz. Therefore we conjecture that similar formula holds true in more general 
situations, in particular for the field theory models. For these models the approach (i) usually 
leads to short distance expansions for the correlation functions, while the approach (ii) gives 
their long distance expansion. The problem of putting into explicit (analytic) correspondence 
these two regimes has been the subject of several works along the last fifteen years (see e.g., 
|4Ul 1411 1421 ESI)) remaining however up to now an open problem. We hope that the method 
presented here could ultimately shed some new light on these topics. 

The paper is organized as follows. In the next section we give the main notations and 
definitions of the XXZ model using the framework of the algebraic Bethe ansatz. In Section 3 
the twisted transfer matrix 7^ is introduced. We describe the properties of this operator and 
related objects, like twisted Bethe equations, their solutions and the scalar product formulas for 
the eigenstates of 7^ with arbitrary vectors. In Section 4 we derive the master equation. The 
last two sections are devoted to the above mentioned two evaluations of the multiple contour 
integral defining this master equation. In Section 5 we reproduce from it the multiple integral 
representation for the a^ correlation function obtained in [^. This representation is given for 
the finite chain and in the thermodynamic limit both for the massless and massive regimes. 
In Section 6 we obtain the form factor type expansion for the a^ correlation function. In the 
Conclusion we give several comments on the key points of our methods. We also discuss several 
perspectives of the obtained result and announce the multiple integral representation for the 
dynamical correlation function of the third components of spin, which in fact is the subject of 
a forthcoming publication |49j . 

2 Algebraic Bethe ansatz for the XXZ chain 

In the framework of the algebraic Bethe ansatz the Hamiltonian (|1.3|) can be obtained from the 
monodromy matrix T(A), that in its turn is completely defined by the i?-matrix. The i?-matrix 
of the XXZ chain acts in the space C^ <^ C^ and is equal to^ 
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sinh(A + rj) J 



coshTy = A. (2.1) 



It is a solution of the Yang-Baxter equation. Identifying one of the two vector spaces of the 
i?-matrix with the quantum space Tim, we obtain the quantum L-operator at the site m 

L^(A) = i?om(A-r//2). (2.2) 



'^Note that we use here a different normalization for the i?- matrix compared to |.-il| . 



Here Rom acts in C^ ® 'Hm- Then the nionodromy matrix T(A) is constructed as an ordered 
product of the L-operators with respect to all the sites of the chain 

^W=(c!aJ o;^|)=iM(A)...L.(A)L.(A). (2.3) 

The operator T(A) acts in V ®Ti, where y -^ C^ is usually called auxiliary space of T(A). The 
Hamiltonian (|1..S)) at /i = is related to T(A) by a 'trace identity' 



dX 



+ const. (2.4) 



A— "J 



Here 

T(A) = tr r(A) = ^(A) + D{\). (2.5) 

Later on we shall consider the inhomogeneous XXZ model, for which 
Lm{\) = Lm{\U) = R^m{\-U), T{X) = Lm{\, ^j) . . . L2{X, ^2)Ll{X, ^l) , (2.6) 

where ^m ai'e arbitrary complex numbers attached to each lattice site that are called inhomo- 
geneity parameters. In the homogeneous limit ^m = ??/2 and we come back to the original 
model. 

The commutation relations between the entries of the monodromy matrix are defined by the 
Yang-Baxter quadratic relations, 

Ri2{Xi - A2)Ti(Ai)r2(A2) = T2(A2)Ti(Ai)i?i2(Ai - A2). (2.7) 

The equation (|2.7j) holds in the space Vi (S>V2 ^TC (where V^- ~ C^). The matrix Tj{X) acts in 
a nontrivial way in the space Vj ^ 7i, while the i2-matrix is nontrivial in Vi V2- 

In the framework of the algebraic Bethe ansatz an arbitrary quantum state can be obtained 
from the states generated by the action of the operators B{X) on the reference state |0) with 
all spins up, 

AT 

\i;) = l[B{X,)\0), N = 0,l,...,M. (2.8) 

i=i 
The eigenstates of the transfer matrix T(/n) can be constructed in the form (|2.8j) . where the 
parameters Xj satisfy the system of Bethe equations 

N N 

a{Xj) Yl sinh(Afe - A^- + 77) = d{Xj) JJ sinh(Afe - A^ - r?), j = l,...,N, (2.9) 

fe=i fc=i 



and a(A), d{X) are the eigenvalues of the operators ^(A) and -D(A) on the reference state. In 
the normahzation (|2.1|1 - (|2.2)) . we have 



M 



a(A) = J]sinh(A-e„ + r/), (2.10) 

a=l 

M 

d{X) = l[smh{X-Ca). (2.11) 

a=l 

Remark. Generahy speaking the system of equations (|2.9|) is neither a necessary nor a sufficient 
condition for the vector (|2.8)) to be an eigenstate of the transfer matrix. For example, the state 
(|2.8j) with all spins down (A^ = M) is the eigenstate of T for generic {A}. On the other hand the 
system (|2.9j) possesses solutions, which do not correspond to any eigenstate of T. We discuss 
all these questions in more details in Section |31 and Appendix 1X1 

The eigenvalue r(/z|{A}) of the operator T{fi) corresponding to an eigenstate of the form 
dUHl) is 

^ \j\\\ ( N TT !i^^M:?^^~ii±^ _L ^^ -pi" sinh(^i-Afc + ry) 
^J-^ smh(Afc - /i) ^^_^ smh(/i - A^) 

The dual states can be constructed similarly to ()2.8() via the operators C(A) 

AT 

(V'l = (0|n^(^i)' iV = 0,l,...,M. (2.13) 

i=i 

Here (0| = |0)+ and the dual eigenstates of T(/i) are given in the form ()2.13|) . where the 
parameters \j satisfy the same system of equations (|2.9() . Generically {4'\ / IV')^ due to the 
involution C(A) = iti?' (A). If, however, the state (|2.8|) is the ground state of the Hamiltonian 
()1.3|) . then {il)\ = c |^)T with c = 1 for A > 1 and c = —1 for —1 < A < 1 respectively. 

In the end of this section we give the explicit representations of the local spin operators in 
terms of the entries of the monodromy matrix. Such a representation is given by the solution 
of the quantum inverse scattering problem j25| I27j 

< = n '^(^^) • Hni,)'y^) ■ n '^-'fe)- (2-14) 

k=l k=l 

Here a" in the r.h.s. acts in the auxiliary space of T'(A), while cr" in the l.h.s. acts in the 
local quantum space Tij. The formulas of the quantum inverse scattering problem permit us to 
embed the problem of calculation of the correlation functions of the local spin operators into 
the algebraic Bethe ansatz. 



3 Twisted transfer matrix and scalar products 

In this section we introduce new objects and notations, which will be used through all the paper. 
Starting from this section we consider only the subspace tyv™/^-^) of the quantum space 7i 
with fixed (but arbitrary) number of spins down N . 

First, we introduce the 'twisted transfer matrix' ^^(/u) 

%{^l) = A{^l) + KD{^x), (3.1) 

where k is a complex parameter. We denote also as IV'sd-^})) a-iid (V^redA})! the eigenstates 
(respectively the dual eigenstates) of the operator '7^(;u) in the subspace JiK^/'^-N) _ The cor- 
responding eigenvalue is denoted as Tk(^|{A}). The same notations without a subscript k 
correspond to the case k = 1 considered in the previous section. 
The eigenstates of Xc(/^) (and their dual states) have the form 

N N 

\M{m = n^(^^-)lO), (^.({A})| = (0|n^(Ai), (3.2) 

with an eigenvalue 

{ \!\r\ ( X TT sinh(Afc - /i + 77) -p-r sinh(^ - Afc + r/) 

^«(^|{A}) = a(^) M — T-rTT ^ + Kd(//) M — —- -^ , (3.3) 

^^_^ smh(Afe - n) ^J-^ smh(/x - Afc) 

where the parameters {A} satisfy the system of twisted Bethe equations 

N N 

aiXj)Y]_s^^H^k- >^j + v) = f^d{Xj)Y\_s'mh{Xk- Xj -r]), j = l,...,N. (3.4) 

It is also convenient to introduce the function 

TV 

y^{f^\{X}) = n sinh(Afc - fi) ■ T«(^|{A}) 

k=l 

N N 

= a{ij) 1 [ sinh(Afc — ^ + r]) + k d{fi) 1 T sinh(Afe — fi — rj) . (3.5) 

fc=l k=l 

In terms of this function the system (|3.4() reads 

y.{X,\W) = 0, j = l,...,N. (3.6) 

Just like in the case k = 1, for generic k not all the solutions of the system (|3.6|) correspond 
to eigenvectors of Tf^ip)- A brief sketch of the properties of the solutions of H3.6|) is given in 
Appendix (for details we refer the reader to the original paper jH]). Here we merely recall 
some definitions and one of the main results of jll] . 



Definition 3.1. A solution {A} of the system iS. 6\) is called admissible, if 



N 



d{X,) J] smh(Aj - Afc + 77) / 0, j = l,...,N, (3.7) 



k=l 



and unadmissible otherwise. A solution is called off-diagonal if the parameters {A} are pair-wise 
distinct and diagonal otherwise. A solution is called degenerated, if the Jacobian of the system 
/ 1 5'. 61) vanishes, and non- degenerated otherwise. A solution is called trivial, if the corresponding 
state IV'red-^})) ^s the zero vector. 

One of the main result of |11] is that, for generic k and {^}, the eigenstates of Xt(/u) corre- 
sponding to the admissible off-diagonal solutions of (|3.6jl form a basis in the space 7-^(W2-a^)_ 
Generically this basis is not normalized and non-orthogonal, although it is orthogonal to the dual 
basis. We would like to mention that for some specific choice of k and {^} certain unadmissible 
solutions may also contribute to the basis of the eigenstates in '}{{m/2~N) ^ 

Consider now the scalar products of eigenstates (|3.2|) and arbitrary states of the form (|2.8j) . 
The explicit results for such scalar products at k = 1 were obtained in |45l I25j . Applying the 
methods used in these papers one can easily prove 

Proposition 3.1. Let {A} satisfy the system ici.b}) . {/u} be generic complex numbers. Then 

N N 

N 



n diXa) ■ X^Hili}, {A}) • d^t f _r,(^,|{A}) j . (3.8) 

a=l \ J / 



Here A'^vd/u}, {A}) is the Cauchy determinant composed of the parameters {A} and {^} 

N 

n sinh(Aa - \h) sinh(/ife - /i^) 

^.(W. (A}) ^ det„ (-sj;;^) - '" . ■ (^■») 

11 smh(//fc-Aa) 

a,b=l 

To make these formulas more explicit we introduce for arbitrary positive integers n and 
n' (n < n') and arbitrary sets of variables Ai,...,A„, /ii,...,//„ and i/i,...,f„/, such that 
{A} C {v}, the following n x n matrix Qf^{{X}, {/i}|{z^}) 



a=l 



{0.^)jk{{X},{fJ-}\{iy}) = a{fik) t{Xj , fik) JJ sinh(z^a - ^fc + ??) 

n' 

Kd{nk)t{nk,Xj) JJsinh(i/a -/Ufe -r/), (3.10) 



a=l 



with 

sinh 77 

*^'^ =— T7^ T— TTT 7- 3.11 

^ '^' sinh(A - ii) smh(A - ^ + v) 

Then the equation f^'A.H^ reads 

N N 

j=i j=i 



rAf 



n;=l d{Xa) 



•detfi,({A},M|{A}). (3.12) 



TV AT 

n sinh(Aa - Afc) sinh(/xb - fia) 

a>b 

Taking the limit {fi} — > {A} in (|3.12)) we obtain the 'square of the norm' of the eigenstate (recall 
that generically {MW)\ + lV'K({A}))t) gagSj 

(^.({A})|V'«({A})) = ^^"=^'^^^"^ • d^t 17«({A}, {A}|{A}). (3.13) 

n sinh(Aa - \h) 

a, 6=1 

This equation can also be written in terms of a Jacobian using 

d^tO.({A},{A}|{A})=det(-^3^«(A,|{A})), (3.14) 

where {A} satisfy the system (|3.()j) . 

4 Master equation for g'^ correlation function 

In this paper we consider the generating functional for the correlation function of the third 
components of spins (o'fo"^_|_]^). Following the papers j46| I39j we define for any complex number 
K the operator Ql^ as^ 

Qi,"^=n(^+^-<)- (4-1) 

n=l ^ ' 

The generating functional is equal to the expectation value 

^^1'™^ (V'({A})|V({A})) ' ^^-'^ 

where |V'({A})) is an eigenstate of T^ii) in the subspace 7-^(*^/2-^). Taking the second 'lattice 

derivative' of {Qim) ^^^ the second derivative with respect to k at k = 1 we extract the 

^Setting K — e" one has Ql^rn = exp (^ X]™=i(l ^ '^n)) and, hence, this is exactly the operator considered in 

ED 



two-point correlation function of the third components of local spins: 

(4.3) 



1 d^ 

■((1 — (yi)[\ — Cr^+i)) = ;^~2 U^l,m+l ~ Ql,m ~ Q2,m+1 + Q2,m)) 



2'' ''' ""+'"' dn'^ 



K = l 



The main result of this paper is an integral representation for the generating functional 

(Ha). 

Theorem 4.1. Let the inhomogeneities {^} be generic and the set {A} be an admissible off- 
diagonal solution of the system \2.y\) . Then there exists kq > such, that for |k| < kq the 
expectation value of the operator Ql^ in the inhomogeneous finite XX Z chain is given by the 
multiple contour integral 



N , N 



r.r^i,,T.rM 7 = 1 a,b=l a=l ^^"-^^ ^' 



r{€}ur{A} ^-1 



Jv 

n 

a=l 



dzk 

n3^.(^a|{4)-det^(^^%P 



■ (4.4) 



The integration contour is such that the only singularities of the integrand \4-4\) within T{^} U 
r{A} which contribute to the integral are the points {^} and {A}. 

We call (|4.4() the master equation. 

Remark 1. The master equation H4.4() gives the expectation value (Qi^) with respect 
to an arbitrary eigenstate iV'd^})) of '^- Iii particular one can chose {A} such that in the 
homogeneous limit ^a -^ ^/2 the corresponding eigenstate lipii^})) goes to the ground state of 
the XXZ Hamiltonian. 

Remark 2. The master equation (|4.4|) is an integral representation of the expectation value 
{Q1m)i which holds true at least for |k| small enough. On the other hand this expectation value 
is evidently a polynomial in n of degree m. Therefore the representation (|4.4j) can be easily 
continued in k from any vicinity of the origin to the whole complex plane. This does not mean, 
however, that one can set k to be an arbitrary specific value directly in the integrand of (|4.4j) . 

Proof oi Theorem \4.1\ In order to compute the expectation value of Ql^, we express this 
operator in terms of the twisted transfer matrix. Due to ()2.14|) one has 

m m 

QU = J{%{ia)X{T-\i,). (4.5) 

a=l b=l 

Then we can use the explicit formula for the multiple action of HoLi '^nixa) for an arbitrary set 



10 



N 



of complex numbers {x} on an arbitrary state (0| n7=i C'(Aj), [S 



N m 



j=l a=l 

min (m,N) 

= E E i?;;(K+},{x^_},{A,+ },{A«_})(0| n C'(xa) n C'(^^)' (4-6) 

n=0 {->>}={Aa^}u{Aa_} ^67+ 6ea_ 

{a:} = {a:T,^}U{a;7_} 

a+l = l7+l=n 

where the coefficient R'^{{x-y_^ }, {x^_ }, {Aa_|_ }, {Aq_ }) is given by 

Rn = \ n sinh(Ab - Aa) J| sinh(x6 - Xa) J| J| sinh(Afe - A^) > 

^ a>b a<b aGOi-^ 6GCK_ '^ 

a,6£ci_|_ a,6G-7_|_ 

X n ^''(^«l{^7+}U{Aa_})-detfi«({xT,+ },{A„+}|{x^^}U{A„„}). (4.7) 

ae7- 

Here O^dx-y^}, {Aa_|_}|{x^^}U {Aq_}) is given by expressions p.lUf) in which the sets {A}, {fi}, 
{u} have to be replaced by {x-y^}, {Aq^} and {x-y_^} U {Xa-} respectively. 

The summation in H4.6() is taken with respect to partitions of the sets {A} and {x} into 
disjoint subsets {A} = {Aq,^} U {Aa_} and {x} = {x^,,,} U {x^ }, such that the number of 
elements in the subsets {Aq,,^} and {x^^} coincides and is equal to n. 

In the paper [HJ we directly applied the equation (|4.6|1 for the computation of (Qi^)' 
specifying Xj = S,j for j = 1, . . . ,m. The peculiarity of the inhomogeneity parameters is that 
d{^j) = 0, what simplifies the formulas. In the present paper, however, we are going to keep 
the parameters {x} arbitrary, i.e. we shall consider the expectation value of the operator 

m m 

Q^{k, {x}) = W %{xa) ■ n '^'H^b), (4.8) 

a=l b=l 

where the parameters xi, . . . , Xm are generic complex numbers. 

The normalized expectation value of Qm(^){x}) on the finite lattice with respect to an 
arbitrary eigenstate of the transfer matrix |^({A})) has the form 

m m 

(M.{mnT,w)-UT-HxMi{\})) 

'^""'- '^»' ^ — "'(^({A))i:^A))) ■ (") 

The action to the right of the product of T^^(x;,) produces merely a numerical factor: 

m m 

'[\T-\x,)\^j{{X])) = X\t-\x,\{X]) . |^({A})). (4.10) 

6=1 b=l 

11 



Acting to the left with the product of T^{xa) by means of (|4.6|) with m < N, and using the 
expression (|l-i.l2j) of the scalar product at k = 1, we obtain 

m 
{Qm{K,{x}))=J2 E A:{{x^^},{x^J,{X^^},{X^_}), (4.11) 

n=0 {a}={Aq_,_}u{Aq_} 

|a+l = l7+l=" 

with 

m -. 

K = \{{r{xb\{X]) JJ sinh(Aa - Xb) W siuh{Xa " Xfo)} J| y^{Xa\{x^+] U {K_}) 



&=1 a&ct— 1S7+ a£7_ 



detjvll({A},{x^^}U{Aa_}|{A}) 



X detS2K({x^+},{Aa_^|i{2;^_^}U{Aa_|) • . ^ ^-^m r\ i , , r\ TTTTvT- ("^-l^) 



Here 3^ and Vt without subscript k are equal respectively to y^ and $7^ at k = 1. The elements 
in the sets {x^^} U {Aq,_} and {Aq,^} U {Aq_} are ordered accordingly. 

Like in pj, we can now perform a re-summation over the partitions of the set {x} by 
introducing contour integrals over auxiliary variables zi, . . . , z„: 



1 / TT dZi 



(Qn.(-,w))=E E n n smh-HA,-x.).^ y ns 

1 fr3^-(^alWu{A„_}) 



n=0 {A} = {Aa_,_}U{Aa_} a=l bSO- PI 1 j = l 



n m ^ m 

V 

■ ■ ■ L sinh(Za - j;i,) ^J^ T{Xa\{X}) ^J^yK{Za\{z}U{Xa,}) 



TT TT J^ TT J^KV^all^J '-^ l^g-j; TT 

a=lb: 

^^.oai/x ii/iM/i 1, detjvl7({A},{4u{A^,}|{A}) 

X detf7.({.},{A.,}|{.} U {A._}) . det^^(|;,|^|;,^^}u{A._}|{A}) ' ('•^') 

where the closed contour T{x} surrounds the points^ xi, . . . ,Xm and do not contain any other 
pole of the integrand, i.e. zeros in C" of the functions 3^K(2;a|{-2} U {Aa_}). Since at this 
stage of the computations xi, . . . , Xm are generic complex numbers, we can always choose them 
separated from the zeros of yK{za\{z} U {Aq_}) for any subset {Ao_}. 

Thus, the sum over partitions of the set {x} in ()4.1ip is replaced with a contour integral. 
In paper j^ we replaced in the thermodynamical limit the sum over partitions of the set {A} 
with the integrals over the support of the spectral density of the ground state. Here we treat 
this sum in a different way. Namely, we perform a second re-summation over the partitions of 

^More precisely, r{a;} is the boundary of a set of polydisks T>a(r) in C". Namely, r{a;} = U^i7?a(r), where 
'£'a{r) = {2 e C" : \zk — Xa\ = r, k = 1,. . . ,n}. The integration contour r{^} U FjA} in 14.411 should be 
understood in a similar manner in C^. 
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the set {A} in (|4.13j) also in terms of a contour integral. Indeed, one can easily see that 



Res 

{z„ + l,...,Zjsl} = {Xa_} 



det ^fc{{zi, ..., Zn}, {Xa+ } U {Xa-}\{zi,. . . , zn}) 



= TT yK{Xa\{zi,---,Zn}U{Xa_})-detn^{{zi,...,Zn},{Xa+}\{zi,---,Zn}U{Xa_}). (4.14) 
a€a_ 

This enables us to express the generating functional {Qm{n, {x})) for the finite XXZ chain as 
a single multiple integral: 



where the closed contour T{x} U r{A} surrounds the points xi, . . . , Xm and Ai, . . . , Xn (spectral 
parameters corresponding to given eigenstate of T{iJ,)) and no other poles of the integrand. 

Finally, in order to reproduce the expectation value {Qi „ , ) we should set the generic 
parameters {x} in (|4.15j) to be equal to the inhomogeneities ^i, . . . ,^m- This definitely can be 
done if the system 

y^iz,\{z}) = 0, j = l,...,N (4.16) 

has no solution inside the integration contour r{^} U TjA}. Thus, we need to analyze the prop- 
erties of the solutions of the system of the twisted Bethe equations. Hereby, since {Qm[f^i {x})) 
is a polynomial in k, it is sufficient to determine this polynomial in an open ball around some 
specific value, for example for |k| small enough. 

The detail analysis of the system 1)4. 16() was done in the paper ^]. We give some of the 
basic statements of this paper in Appendix \M Here we present only the results necessary for 
writing and evaluating the integral (|4.15j) . 

Lemma 4.1. Let the inhomogeneities ^i,---,^m ^e generic, Xi,...,X]\f be an admissible off- 
diagonal solution of \2.fJ\) and \k\ be small enough. Then all admissible off-diagonal solutions 
of the system J^. J6] j are separated from the points {^} and {A}. 

Corollary 4.1. There exists a contour r{^} U r{A} such that the points {^} and {A} are 
inside this contour, while all admissible off-diagonal solutions of the system i4-l^\ ) '^'"^ outside 

r{e}ur{A}. 

Lemma 4.2. Let the contour T{^} U r{A} satisfy the conditions of the corollary \4 1[ and 
Xk ^ ^k- Then: 

1) the only poles inside the contour T{£f\ U FjA} which provide non-vanishing contribution 
to the integral | |^.i5| ) are in {^} U {A}; 

2) the only poles outside the contour r{^} U r{A} which provide non-vanishing contribution 
to the integral | |^.i5| ) are the admissible off-diagonal solutions of the system | |^.i61 ). 
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The proof of these lemmas is given in the Appendix fusing the results of |44j . 

Lemma 14.21 guarantees that for Xk = S,k there are no any other poles of the integrand (|4.15|) 
inside the contour r{^} U r{A} and contributing to the integral except Zj = ^^ and Zj = Afc. 
Thus, setting x^ = ^k and using d{S^k) = we arrive at 



N , N m 



Wi.m/ ^, jP ii 27ri ■ 11 11 V sinh(z„ - ^,) ' sinh(A, - 6 + r/) ^ 

detiv ^n{{z}, {X}\{z}) • detiv ^({A}, {z}\{\}) 



r{5}ur{A} ^-^ °—^ ^-^ 



N 

Uy^iza\{z})-detNni{X},{X}\{X}) 

a=l 



(4.17) 



It remains to express the determinants of the matrices Q^ and fi in terms of Jacobians using 
(|3.8() - (|3.14() . and we obtain the master equation H4.4() .D 

5 Multiple integral representation 

Lemma 14.21 permits us to evaluate the integral 1)4.41) in two different ways: either to compute 
the residues of the integrand inside the contour r{^}ur{A}, or to compute the residues outside 
this contour. In this section we consider the first way, which immediately leads us to the 
representation for (Ql^) obtained in |3T1. 

For this purpose it is more convenient to use (|4.17j) . which in fact is equivalent to (|4.4j) . 
This multiple integral can be presented as 

N N ^ n ^ N-n 






r{€}ur{A} •'-^ "-^ r{c} ■^'-^ r{A} ■^'-^ 

Hereby, since the number of the poles inside r{,^} is m and the integrand vanishes as soon as 
Zj = Zk, the sum in (|5.1)) is actually restricted to n < m. Evaluating N — n integrals in the 
points {Xa_} we arrive at ()4.13() with Xk = S,k, and hence d{xk) = 0. Moreover, we can also set 
d{zk) = 0, since the remaining integrals surround only the poles at Zk = S,j and after evaluation 
of these integrals the functions d{zk) vanish. Finally, using the fact that the set {A} satisfies 
the system 1)2. 9|l we obtain 

, , .., n c(Ai,i{4, {>.-)■{?}) 

n=0{A}={A,+ }u{A,_} -plf) J=l UUizh\{z},{Xa_},{^}) 

l"+N'i ^^^ 6=1 

det^f^({A},{4u{A„_}|{A})| 

xdetiV4{A»J,{4 o(JMn InSX 11/^ ' ^'^ 

n + detAriZ(|A},{AQ^} U|Aa_|l|A}) 
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where 



sinh(z^ - ^a) 



U{u\{z},{Xa_},{0) = a{u)T\smh{za-u+v) n sinh(Aa-z/+r/)rr . ^ v 

-'■-'■ -'■-'- -'■-'- smniz^ — to + 77 

a=l aGa- a=l ^ sa i /; 



and 



(MK)ifc({A}, {z}) = t{zk, \j) + Kt(Aj, Zfc) J]^ 



sinh(Aa - Xj + i]) sinh(Aj - Za + rj) 



L sinh(Aj - Xa + v) smh{za - Xj + rj) ' 



(5.3) 



(5.4) 



Observe that, unhke in the master equation, the integrand in the r.h.s. of (|5.2|) is a polynomial 
in K, since n enters only the determinant of the matrix M^ (|5.4|) . Therefore one can set k to be 
an arbitrary complex number directly in the integral representation (|5.2j) . 

This representation also allows one to proceed, if necessary, to the homogeneous limit simply 
by setting ^k -^ ilf^- Finally, the equation (|5.2|) is convenient for taking the thermodynamic limit 
A^, M — > oo N/M = const. We refer the reader for the details to the paper [HJ and here simply 
recall the final result for the ground state expectation value of {Qim) ^^ ^^^ thermodynamic 
limit and in an external magnetic field in the massive and massless regimes: 



(Qlm) = E 



ra=0 
n n 

nn 



{niy 



dz^ 



r{?} 



i=i 



nsy/"^nn 

„_1 J Ly L_i „_i 



=la=l 



sinh(Ab - ia) sinh(zfe - ja + ri) 
sm\i{zb - ia) sinh(Ab -ia + v) 



sinh(Aa - Zb + T]) sinh(zj, - Ag + 7?) 
^Jl ^j^ sinh(Aa - Xb + r]) sinh(za - Zb + r]) 



■ detr 



M.i{X}\{z}) 



detr. 



p{Xj,Zk) . (5.5) 



Here the function p{X, z) is the 'inhomogeneous spectral density' of the ground state, having 
the support on the contour C as defined in pT] . 

6 Form factor expansion 

Evaluating the integral (|4.4() by the residues outside the integration contour we arrive at the 
expansion over form factors for {Ql^)- Recall that due to Lemma EH the only poles outside 
the contour r{^} U r{A} which contribute to the integral (|4.4|) . are the admissible off-diagonal 
solutions of (14.161). Hence 



N 



\Ql,m) 



{fi} a, 6=1 



MfcJ 






n 



det 



TV 



dfJ-k 



detjv ('M^^l 






detJ 'y±\^^^A detJ'-^^W^ 



dfi 



V 



(6.1) 



where the sum is taken with respect to all admissible off-diagonal solutions fii, . . . , pN of the sys- 
tem ()4.16() . Due to the Theorem lA. li the Jacobian matrix ^^^^(^fcll^D/t^Mj is non-degenerated. 
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Using the formula (|3.8j) we identify each Jacobian in (|6.1|) with the corresponding scalar 
products, leading to 

.^. ,_vft T.felM) (V'({A})IV'.(M)) {^.{{i^})m{\})) ._. 

It remains to use that the state \ipK,{{lA)) is the eigenstate of X<;(0 with the eigenvalue Tk(.^|{/u}) 
and the state IV'd^})) is the eigenstate of T(^) with the eigenvalue t(,^|{A}). Thus, the equation 
(|6.2j) can be written in the form 

(V'dA})! n ^.(6)i^.(M)) • (^«(M)i n ^-HCb)iV'({A})) 

^^^•™^ ^ ^ (^k(M)IV'.(W))-W{A})1V'({A})) • ^^-^^ 

Observe that we did not use the statement b) of the Theorem lA.ll on the completeness of the set 
IV^Kd/i}))- The sum over eigenstates of 7^ appears automatically as the result of the evaluation 
of the multiple integral (|4.4)) by the residues outside the integration contour. 

Taking the second lattice derivative of (jH^ and then differentiating twice with respect to k 
at K = 1 (see AppendixEl) we obtain the form factor type expansion directly for the correlation 
function of the third components of the spin 

/ - ^ \-/ ^N / ^ \^ \- (^({A})ki1V'(M)) • (^(M)K+il^({A})) ,„ .^ 

(-1-^+1) - (-1) • (-^+1) + 1^ WM)|V'(W))-(^({A})|^({A})) ' ^'-'^ 

where the form factors of a'^ are given by ()B.lfl|) . Observe that here the summation is taken with 
respect to the eigenstates of the operator T, but not 7^. In other words, in the homogeneous 
limit, this is the sum over the excited states of the Hamiltonian. 

Of course, the equation (|6.4|) might be obtained by means of inserting the complete set 
of the eigenstates of the Hamiltonian between local spin operators. Such a representation, 
however, would be quite formal without a precise description of the set of the states on which 
the summation has to be performed. On the one hand it is clear that the sum should be taken 
with respect to solutions of the system ()2.9() . but on the other hand it is also clear that not 
any solution of this system corresponds to an eigenstate of the Hamiltonian. Note that one 
can not simply restrict the sum in ()6.4() to the admissible off-diagonal solutions, since in the 
homogeneous limit certain unadmissible solutions of the system (|2.9j) give rise to the basis of 
the eigenstates. 

The master equation approach gives the way to overcome these difficulties. Namely, if we 
use the approach (|1.2j) . we can insert the set of the eigenstates of the twisted transfer matrix 
between local operators. Then the summation in (|6.4j) is restricted to the admissible off-diagonal 
solutions of the system (|3.6() . The unadmissible solutions corresponding to the eigenstates of 
the Hamiltonian appear in this case only as the limit of admissible solutions at k ^ 1. Since for 
|ac| small enough all admissible solutions are in the vicinities of the points {^ — 77}, the sum over 
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excited states can be written as a multiple integral with respect to a contour surrounding these 
points. This allows one to obtain the master type equation directly for the spin-spin correlation 
functions, starting form their form factor expansions. 

7 Conclusion 

The main result of this paper is the master equation (|4.4j) . We have shown that this equa- 
tion draws a link between multiple integral representations and form factor expansions for the 
correlation functions. 

The master equation sheds a new light on the role of the auxiliary contour integrals in the 
representations for the correlation functions derived in |HJ. Originally they appeared as a result 
of a re-summation of the elementary blocks obtained in j22l 1231 I26j . In the framework of the 
master equation approach the same auxiliary contour integrals are equivalent to the sum over 
excited states. 

Using the master equation method one can derive multiple integral representations for other 
correlation functions of the XX Z model, including many-point correlators. Indeed, an arbitrary 
correlation function in the finite chain can be reduced to the multiple sum over complete set of 
the eigenstates of the twisted transfer matrix. On the other hand the explicit formulas for the 
form factors of the local spin operators obtained in [2S1 can be easily generalized for the case 
when one or both states are eigenstates of the twisted transfer matrix. Presenting this sum as a 
contour integral of the type (|4.4() and evaluating this integral by the residues inside this contour 
we can express an arbitrary correlation function of the XXZ model as a multiple integral of 
the form (|5.5|) . 

Moreover, we conjecture that representations similar to 1)4. 4|) should exist for the correlation 
functions of other models solvable by the algebraic Bethe ansatz. This conjecture is based on 
the fact that the integrand of (|4.4j) depends mostly on the eigenvalues of the twisted transfer 
matrix, which is a typical object for the algebraic Bethe ansatz. As we have mentioned already 
in Introduction, it would be very desirable to obtain an analog of this master equation in the 
case of the field theory models. For these models such an integral representation could give an 
analytic link between short distance and long distance expansions of the correlation functions. 

Finally, the master equation method opens a way to obtain multiple integral representations 
for time-dependent correlation functions. This subject will be considered in our forthcoming 
publication |1^. Here we only announce a generalization of the representation 1)5. 5() for the 
time-dependent case in the massive and massless regimes 



K(OX+i(i)) = 2K(0)) -1 + 2DI ^Qn{m,t) 



(7.1) 

K=l 
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where D^ means the second lattice derivative and 



n (rr, +\-\^ _J_ [ rrn\ I TT ^ TT sinh(Aa -Zb + T]) smhjzb - Ag + r?) 
J.[m, t) - 2^ ^^,^2 J ^ ^ f [I 2^, • 11 sinh(A, -Xb + v) smh(z, - z, + r?) 

Tl — U /-# -n (- I 77-1 7 — J- (2,0 — i 



c r{±f } ^ 



n 
X 
6=1 



YT^it{Eiz,yEix,))+imipiz,)-p(x,)) detM«({A}|{z}) •det[7^"(A.,Zfe|{A},{z})]. (7.2) 



Here M^ is given by (|5.4I) . The functions i?(A) and p(A) are the bare one-particle energy and 
momentum 

^r.^ 2sinh2r/ ,,, , /sinh(A-2)\ 

^^^^ = sinh(A+i)sinh(A-i)' ^^^^ = ^ ^°^ Uinh(A + i)J ' ^^'^^ 

The contour T{±7]/2} surrounds the points ±??/2 and does not contain any other singularities 
of the integrand. The function 7^^(A, ^[{A}, {z}), as a function of z and other arguments fixed, 
has a cut between the points z = r]/2 and z = —rj/2. Therefore it is defined differently in the 
vicinities of these points 

{p{\, z), z ~ ?7/2; 

.--lnC\ 7^r,\ n" sinh(z~At,+>;)sinh(2fc-z+r?) /„ , (7.4) 

-h, /5iA, Z + rj)- [ 1^^^ sinh(Af,-2+»?) smh(2-26+»?) ' Z^-rj/Z. 

where p(A, z) is the inhomogeneous spectral density of the ground state. 

Observe that the time and the distance dependence of the generating function (|7.2|) is asso- 
ciated to the bare energy -E'(A) and momentum p(A), making the equation (|7.2j) very suggestive. 
In the limit t = the integrand in H7.2|) is a holomorphic function in the vicinity of —ri/2, there- 
fore the integrals over T{—r]/2} vanish, and we arrive at the time-independent representation 
(jESJ. 
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A Solutions of the twisted Bethe equations 

In this Appendix we prove Lemmas 14.11 14.21 using the results of j44j . Consider the system of 
equations (|4.16jl 

y^{z,\{z})={), j = l,...,N. (A.l) 

Without loss of generality we can identify two solutions of this system {z} and {z'} if they are 
equal modulo ivr. 

It was proved in 44^ that 

Theorem A.l. '^44j Let k and inhomogeneities {S,} be generic. Then 

a) All admissible off-diagonal solutions of the system HA . 1\) are non-degenerated. 

b) The set of the states \3.'J\) corresponding to the admissible off-diagonal solutions form a 
basis in the subspace 7-^(W2-A')_ 

The strategy of t44^ was to consider the limit of the system ()A.1|) at k — > and then 
to deform the solutions {z(0)} to the case k 7^ 0. Then the statement a) follows from the 
implicit function theorem. Indeed, the solutions of the system ()A.1|) at k = are evident. In 
particular all admissible solutions have form Zj{Q) = ^p. — rj. One can easily check that the 
Jacobian matrix dyK{zj\{z})/dzk at k = for admissible solutions is a diagonal matrix with 
non vanishing entries. Hence, in this case the solution {z{k)} is a holomorphic deformation 
of the solution {2;(0)}, and therefore for \k\ small enough all admissible solutions are in the 
vicinities of {£, — rj}. On the other hand it is clear that any admissible solution at k 7^ 1 is 
separated from the admissible solutions at k = 1. Thus, admissible solutions are separated from 
the points {^} and {A}. This proves Lemma l4. II 

In order to prove Lemma 14.21 we use the properties of the matrix il^dA}, {/i}|{i^}) (|.S.10|) 
formulated in the following two lemmas. Let the parameters Ai, . . . , A„ and /xi, ...,//„ be generic 
complex numbers, k be arbitrary complex. 

Lemma A.l. Suppose there exist fia, fJ-b C {fi}, such that fia = S,pi l^b = (,p — Vj where ^p is 
one of the inhomogeneity parameters. Then det„ il^dA}, {/x}|{A}) = 0. 

Proof. We have d{fj,a) = o-ifJ-t) = and t(/i6, Xj) = t{Xj,fia)- Thus, the columns {0,fi)ja and 
{QK)jf) are proportional to each other, hence det„r2fj = 0. □ 

Lemma A. 2. Let q^ = e^^ be a root of unity, i.e. r] = inQ/P, where Q < P are positive 
integers. Suppose there exists {Am, . . . Aap} C {A}, such that 

sinh(Aa^^-^ - Xa^ + r/) = 0, with Aap+i = Aai. (A.2) 

Thendetnn^{{X},{n}\{X}) =0. 

Proof. It is easy to see that 

p p 

J^%,AaJ = ^t(A,^,,/x) = 0, (A.3) 

i=i i=i 

19 



since both these sums are ivr-periodical holomorphic functions of /i vanishing at ^ ^ ±00. Hence 
the matrix J7k({A}, {//}, {A}) contains hnearly depended hues and its determinant vanishes. n 

Note that if the set {A} contains several subsets of the type ()A.2|) . then the order of the zero 
of det,i r^fjdA}, {m}|{A}) is not less than the number of such subsets. 

If q^ is not root of unity, then unadmissible solutions of the system HA.1|) contain a pair 
^a = ip and z\, = ^p — rj. Off-diagonal unadmissible solutions are non-degenerated at k = 
|44j . Hence, due to Lemma lA. II thev do not give a contribution to the integrals (|4.17|) . ()4.15() . 
The diagonal solutions {z{k)} of (jA.lj) obtained as a deformation of diagonal solutions {-z(O)} 
at K = preserve their multiplicity jH]. Hence, they also do not give any contribution to 
the integrals (|4.17|) , (|4.15j) due to vanishing of det jv ^k and detjv ^ as soon as two or more z 
coincide. 

If q"^ is a root of unity, then unadmissible solutions may contain 'strings' {za } satisfying 
the condition (|A.2|) . Such solutions are not isolated since the value of Za^ is not fixed. We have 
seen however that due to Lemma lA. 21 the determinant of the matrix ^^^{{z}, {A}|{z}) vanishes 
exactly on these string type solutions. Hence, these unadmissible solutions of the system ()A.1|I 
do not contribute to the integrals (|4.17l) . (|4.15|) . 

Thus, the only solutions of the system (|A.1|) which can give non- vanishing contribution to 
the integrals (|4.17l) . ()4.15|) are admissible off-diagonal solutions, which are in the vicinities of 
{C ~ v} foi^ I'^l small enough. This proves Lemma 14.21 

B The a^ form factor 

The explicit formulas for form factors of local spin operators in the finite XX Z chain were 
obtained in [5S]. Here we propose slightly modified method to derive the form factor of the 
operator o"^. 

Consider a matrix element of the operator Qi^ between an eigenstate iV'dA})) of the 
transfer matrix T(z^) (for example, ground state) and an eigenstate (V'/tdA'})! of the operator 
7^(i^). Using (|4.5|) we immediately obtain 

(MiK'^mQlmmim = n ^^^f}\nx~!^ ■ (V'.({^('^)})!V'({A})). (B.i) 

Here we have indicated explicitly that the parameters {/i} depend on k, for they are solutions 
of the system ()3.6|l . On the other hand it is clear that 



^(v«({^w})i[Q?,„+i - QiraW^am 



= ^(V'({Ml)})l(l-<+i)lV'({A})). (B.2) 
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Thus, we have 
1 



d 



,(*mmi(i-<,Mi{m = j-^n .(j^,,,,. 



T-K{^a\Mn)}) f T^iCm+lllfiin)}) 



,N 



Ua=ld'i^^a{K)) 



N N 

n sinh(Aa - Xb) sinh{fib{K) - ^JLai^)) 

a>b 



detn^{{fiiK)},{X}\{fi{K)}) 



(B.3) 



K=l 



In order to evaluate exphcitly the derivative over k in (|B.3|) one should distinguish two cases: 
{yu(l)} = {A} and {^(1)} / {A}. In the first case T^{Cm+i\{tJ'ii^)}) —^ r(^m+i|{A}) as k ^ 1, 
therefore 



(V({A})|V({A})) dK\ t{U+i\W) 



K=l 



fc=l 



cLk 



■ t{Xk,Cm+l)- 



(B.4) 
(B.5) 



K = l 



The derivatives dijik{n) / dn can be found from ()3.4p via 



N 



E 

fc=i 



dy^i^'MlA) dfj.k{K) 



dfj,k 



dn 



N 



(B.6) 



+ d{^i.j) JJ sinh(^a - ^ij -r]) =0. 

'^=1 a=l 

In the second case we can compute explicitly the derivative of detri^- Indeed, consider the 
A'^-dimensional vector-column v with the components 

N N 

^fc = JJ sinh(;Ufc(K;) - Xa) JJ sinh"^(/Xfc(K) - Ha{n))- (B.7) 



a=l 






If {//(I)} ^ {A}, then this vector at k = 1 has at least one non-zero component, say vjy ^ 0. 
Then multiplying the /c-th column of (J7k)jA: by Vk/vN and adding the first {N — 1) columns to 
the last one, we obtain 

N-l , , N 






l-K 

VN 



i(Aj) JJsinh(Aa-Aj+7?). (Bi 



fc=l a=l 

(see Appendix B of 31^ for the proof). Thus, the last column is proportional to 1 — «;, hence, 
taking the derivative over k of det il^ one has to differentiate only this column. This gives 

d 



^^d^tfi,({M'^)},{A}|{/i(K)}) 



detf^({/x(l)},{A}|{Ml)}), 



K = l 



N 



(B.9) 



where 



(f^),fe(M,{A}|M) 



A; = l,...,A^-l, 



f (i^.U(M,{A}|M)|^^,, 

- iz ■ «(A,) \fa=i smh(Aa - A, + r?), k = N. 
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Thus, for {fi} 7^ {A} we obtain 



(^(M)k..il^({A})) = 2 n ^(^^^ • (1 - ,(e™,,|{A})J 



n-i^^^'^) •detO(M,{A}|M). (B.IO) 



n smh(Aa - Afo) smh(^b - fia) 

a>b 
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